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PREFACE 


Thta document reports on the results of an experimental study 
conducted to determine the geometric emd radiometric effects of 
double resampling ('’bi-resampling”) performed on digital images in the 
process of performing map projection transformations. 
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^Section 1 Introduction nnd siuran-iry 


Thin in th<? finnl r**port for the Experimental Bi-resampled 
Data r.turly nerformed by lii*! for the Doddard Srace Fliftht -'enter under 
contract 'JAH5-J?37 0>3 . 
l.i Overview 

The purnose of this study is to determine the radiometric 
and geometric effects caused bv the double resampling of an image. A 
double resannled { bi-resampled ) image in this experiment is an image 
created bv first resamnling raw, uncorrected data into a Universal 
Transverse ’’ercator (UT'4) nrojection, and then resamnling this UTM 
image into e Lambert Conformal Conic projectio; . This image is then 
subjected to various comparative analyses with respect to singly 
rerampled imagery which was created by resampling the raw data directly 
into the Lambert projection. The two scenes considered in this study are 
l) band 5 of scene E-1080-15192 ( a Chesaperke scene ), and 
C) lands ’» .rough 7 of a sub image of scene E-2183-16433 
covering a LACIE inter, ^ve training site. 

These scenes are designated as scene A and scene B, respectively, 
throughout this report. 

•"he accuracy of uhe manpirig functions used in the resampling 
nrocrss, ro»* both the singly and doubly resampled imagery, is discussed 
i*. Tection P.. fection 3 presents the results obtained in the measurements 
of the registration accuracy between the two versions of scene A. 

Various puantlt.'itive measures of the radiometric degradation due 
to the bl-re.»amnling are discussed in Sectior U. The final secti"'n 
discuss the results obtained for classification processing of 
the tw(. rsions of scene 3. 
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1.2 Summary of Results 


- •■I'l'.' -'i's'].,!: !’! V ■ >!' TIi 
vAL ii; I 

The results obtained in this study sure summarized by the 
following Itens; 

a) The misregistration of the tvo Tersions of scene A was 
foiud to be bounded by 0.29 pixels (the pixel spacing is 50.8 
meters), and for the two ▼ersions of scene B the misregistration 
bound was 0.20 pixels. These are gross misregistration bounds. 

b) IBM's correlation program was applied to 9 widely 
distributed features in scene A, and the maxiimm misregistration 
of any of these features was found to be less ttasm 0.1 pixel. 

The program ms also applied to 19 clustered features. The registration 
results obtained show randcan misregistration errors. 

c) Radioiretric differences on a pixel-by-pixel basis were 
found to be small. Difference imagery, radiometrlcally stretched 
for photographic recording, exhibits differences primarily along 
sharp edges, such as land-water lnterf8u:ee and airport runways. 
Considering all bands of both scenes, roughly 97^ of the pixel 
values in the bl-resampled imagery differed by less than 1 count 
from the corresponding values in the singly resampled imagery. 

d) Edges were degraded slightly in the bl-resampllng process. 

The bi-resampled data displayed "overshoot" on edges, similar to 
the Glbb'e phen^tmenon of Fourier series at a discontinuity. No 
noticeable spreading of edges was observed. 

e) '.he classification results obtained for the tvo versions 
of scene 3 were essentlal.ly the same. Classification accuracy. 
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with respact to ground truth data, is InaanBitive to the singly 
or doubly resampled nature of the data classified. The largest 
discrepancy observed in the areal proportion estimates obtained 
from the two data sets was l.l]^. 
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flection r* flapping Functions 

IB*1 chose to resample, or project, the uncorrected scenes and 
the ITTf scenes to the Lambert Conformal Conic (LCC) projection* The 
resnmnling task requires an inverce mapping function which maps pixels in 
the LCC image space to pixels in the Input image space, which in this 
study was either the uncorrected image space or the UTM image space. These 
inverse mapping functions are denoted by C,, which maps the LCC image 
space to the UTM image space, and H, which maps the LCC image space to 
the uncorrected data. The function G is the inverse mapping function for 
the bi-resampling process, and the function H is the mapping function for 
the single resampling process. 

These functions vere developed as a series of transformations 
ana are illustrated in Figure 2-1. The function G is the mathematical 
composition of the following transformations: 

a) a rotation and scale change which transforms pixels in the 
LCn image snace to (x,y) coordinates in LCC space, 

b) TTiA’^iKPT, which maps (x,y) coordinates in LCC space to geodetic 
latitude and longitude coordinates. Details of this transformation and 
its inverse transformation, LAMBERT, are provided in the appendix. 

c) ULL, which maps geodetic latitude and longitude to the UTfi 
coordinates of northing and easting. This transformation was provided 

bv an APTi version of the equations employed by IBM in correcting imagery 
to the fTTT? projection. 

d) Tp, a translation, rotation and scale change which maps UTM 
northing and easting to l^TM image coordinates. 

The manning ^unction H is the composition of the above function G with F, 
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vhere F represents the Inverse mapping function used in the original 
production of the UTM image. For scene 1080-15192, F is a pair of hivariate, 
fifth order nolynomlals. For the LkClE intensive training site, F is a pair 
bivariate, second order polynomials. The UTM image of the LACIE site 
tra.3 not available, and so a region ' ' the uncorrected image, consisting of 
approximately UoOxUoO samples, ms resampled to create an image of 
this area. The imag‘d size vas chosen so that image edge effects vould 
not be encountered in the subsequent resampling to LCC space. The explicit 
mapping function used in the original creation of the UTM image of scene B 
vas not used to create this UTM subimage. Instead, a polynomial, approximation 
to this explicit mapping function was developed by the method of least 
squares over a region containing Hand County and the LACIE site. 

At this point it should be noted that if the functions G auid H 
were used for mapping etll image coordinates in the resampling program, 
there would be no registration error between the singly and doubly 
resaitmled images. However, it is impracticeil to specify G amd H explicitly 
In the resampling program because of their computationail complexity. 

Instead, the procedure used is to develop a pair of bivariate polynomials 
to annroximate these mapping functions, and a regulair grid as am interx>olation 
net over the image. The resampling logic uses the polynomials to map the grid 
uoints and bilinear interpolation to map other image points. This introduces 
geometric errors which can be characterized as either errors due to the 
nolvnomial approximation (E) or errors due to the bilinear approximation 
of the polynomial at off-grid points (Eg). Bounds for these errors in the 
resampling performed in this study are presented in Table 2-1. The errors 
are expressed as a fraction of the 50.8 meter pixel spacing. For scene B, 
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which was small, polyncmlal approximations ware not used. In this case the 
mapping was provided explicitly. In Table 2-1, the notation ''IX^CCUIM)" 
rfifera to the resampling to the LCC Ima^e from the image and "LCC(II)'' 
efers to resampling to the LCC Image from the uncorrected Image. Generally, 
the errors In the table are pessimistic, since the errors are probably 
not additive, and the bounds are only approached near the comers of the 
Image. 
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E( pixels) 

Eg(plxel« ) 

Or id 8i>acing 

ftcene A, ITRi projection 

0.01 

0.06 

100 pixels 

ICC(11T»0 projection 

0.00 

0.00 

200 pixels 

LCC(ll) projection 

0.01 

0.19 

200 pixels 

Total misref^istration bound 

0.02 

0.25 


Maximum Misregistration 

0.27 pixel 


Scene UTM projection 

NA 

0.10 

<.0<f lines 
samples 

LCC(UTM) projection 

NA 

0.00 

c0.< pixels 

LCC(ll) projection 

NA 

0.10 

e.\j^ pixels 

Total misregistration error 

NA 

0.20 


Maximum Misregistration 

0.20 pixel 



Note: LCC(II) - Lambert derived directly from raw data. 
LCC(UTM) - Lambert derived from UTM iinai;e. 


Table 2-1. Registration Error Bo\inds 
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?>otion 3 M«aBur»*<1 fJaptlBtratlon Rrror 


Misregistration was investigated with IBM's control point 
correletion algoritbn. Nine features widely distributed over scene A, and 
19 closely clustered features in the Cantp Hill region were chosen. The 
locations of these features are indicated in Figure 3-1 « and the 9 widely 
distributed features from the singly and doubly resampled images are 
presented in Figures 3-2 and 3-3. In the latter figures, the features have 
been enlarged to a scale of approximately 1:50000. Hhadeprints of the 9 
features are provided in Figure 3-*». 

The correlation program's registration resxilts are given in 
Tables 3-1 and 3-2. In these tables, line and sample misregistration are 
given is units of the 50.8 meter output pixel spacing, and are designated as 
Ax and Ay, respectively. The worst case line and sample misregistration 
was found to be Ax b o. 07 and Ay = 0.07. The measured standard deviation of 
the misregistration in the cluster area suggest that the errors are randcxa, 
assuming the error function to be slowly varying. Included in these 
misregistration errors are the estimation errors inherent in the cort elation 
program. To estimate these errors, the algorithm was applied to the singly 
resampled scene alone. The resxxlting autocorrelation errors, which axe 


purely algorithm estimation errors. 

are: 



Mean Ax 

Mean Ay 

Std Dev Ax 

Std Dev Ay 

Cluster Region 0.002 

(lo features) 

-0.001 

0.01 

0.02 

nistributed Features 0.002 
(9 features) 

-0.003 




These results indicate that the misregistration results in Tables 3-1 and 
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3-<^ are of the same order of magnitude as the correlation estimation errors. 

In both cases, the errors are negligably small. These control point results 
are well within the -^viously developed bounds on misregistration (Section 2). 
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Feature 

Line Misregistration 
Ax (pixels) 

Sample Misregistration 
Ay (pixels) 

1 

0.02 

0.07 

? 

0.03 

-0.03 

3 

0.03 

C.OU 

u 

-0.01 

-0.05 

5 

-o.oij 

-0.05 

6 

0.00 

0.01 

7 

-0.01 

0.00 

« 

0.00 

0.C2 

u 

0.05 

0.01 

mean 

0.008 

0.002 


Table 3-1. Measured Misregistration Errors for Nine Widely Distributed 
Features in Fcene E-1080-15192 
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Cluster 

Line MlBre/rlstratlon 

Sample Mlseegistratlon 

Feature 

Ax (pixels) 

Ay (pixels) 

1 

0.00 

U.Ol 

2 

- 0.06 

0.00 

i 

0.03 

-0.02 

U 

0.05 

0.02 

5 

0.07 

-O.OU 

6 

0.05 

-0.03 

Y 

-0.02 

-0.01 

8 

-0.01 

0.01 

9 

0.05 

0.01 

10 

U.OO 

0.00 

11 

u.cU 

-0.05 

12 

-0.05 

-0.03 

13 

0.01 

0.02 

ll 

-0.03 

-O.OU 

15 

-0.03 

0.02 

16 

0.01 

-0.03 

17 

-0.01 

0.02 

18 

0.03 

0.00 

• 19 

-0.02 

o.o6 

mean 

O.OOC 

-o.oou 

3td dev 

O.OU 

0.03 


3-2, ^teasured 'lisrefristration Krrors for 19 Clustered Features 
in the Camp Hill Region of Scene E-1080-15192. 



This figure is a photograph contained 
in the pocket at the back of this report. 


Figure 3-1. Feature Location, Cluster Region, and High, Medium and 
I.OW T^requency Regions of Scene E-1080-15192 
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This figure is a photograph contained 
in the pocket at the back of this report. 


Figure 1-2. Nine Features from the Singly Resampled Image of 
Scene E-1080-15192. 
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This figure Is a photograph contained 
in the pocket at the hack of this report. 


Figure 3-3. Nine ’='eat ires fron the Doubly Resampled Image of 
Scene E-1080- 5192 
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flection Radicmietrlc Degradation 


Visually, there appears to be little degradation of the 
bl-resanpled Imsige of scene A. The singly resampled image auad the 
doubly resampled image are shown in Figures *•-! and h-2, 
respectively. A visual image of the radiometric differences was 
produced by subtracting the bl-resampled image from the singly 
resampled image, adding a bias, and exaggerating the radiometric 
scale. This difference image is shown in Figure l(-3. Keeping in 
mind the fact that these differences are slight, several comments 
can be made regarding the difference image. First, one can see in 
the difference image the lemd outlines and features with sharp 
edges, such as ' rports. Second, the overall image tends to a 
uniform gray. Indicating that there sure no gross geometric distortions 
in the registration of the two Images. Sharp edges stand out in 
the difference image due to the "overshoot" of the cubic 
convolution algorithm employed in the resampling process. This 
characteristic of the cubic convolution algorithm is discussed 
in more detail later in this section. 

Quantitative measures of the effects of double resampling 
were obtained through histograms of the difference image, histograms 
and statistical evaluation of corresponding subimages from each 
of the resampled Images, and power spectra analyses of these 
corresponding subimages. 

a) A quantitative measure of the radiometric differences is 
provided in the histograms of the difference images and the two 


U-1 



imageB each of scenes A and B, Figures U-U through U-9« In Figure 
U-9, the histogram of the difference image of the interior region 
of scene A, ^1% of all pixels have count differences of at most 
one count. 

Two comments need to he made at this point regarding the 
UTM image of scene A, from which the hi-resanpled image was created. 
First, the cubic convolution algorithm had a round-off error at 
the time the UTM image was produced. The effect of this problem 
is that the average count In the UTM image is 0.5 count lower 
than the average count in the raw, uncollected data. This is 
supported by the histogram (Figure U-6) which shows the average 
count of the difference image to be lU0.5« Since the bias added 
to the difference image was lUo, this Implies that the singly 
resampled image has an average count 0.5 greater than the 
bl-resampled image, the co'^nt loss in the UTM image having been 
carried over to the bl-resampled image. Second, the UTM image did 
not have a sufficiently large border to eliminate "edge problems" 

when resampling to produce the bi-resampled image. The edge 

/ 

problem has the effect of spreading the histogram of the difference 
image. To eliminate the edge problem effect on the histogram, the 
interior region of the difference image was hlstogrammed, and is 
shown in Flgvire k-g. This figure shows the ellmlnatloa of the 
tails of the histogram in Figure U-6. Neither of the above two 
problems apply to scene B. 

b) Three regions of each scene were selected on the basis 
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of their relatively hifth, medium, and low spatial frequency 
characteristics. Areas selected in scene A are boxed in Flfptre 
3-1, and histograms of each eurea for the singly and doubly 
resampled images are provided in Figures U-IO and b-ll. A 
shadeprlnt of band 5 of scene B is given in Figure U-12, and the 
areas selected in this scene are shown in Figure U-13. 

Radiometric averages and standard deviations were computed 
and plotted as a function of area size for the three regions from 
the singly and doubly resampled Images. The various areas used 
in the computation are centered in their respective regions. Plots 
of these for scene A and the four bands of scene B are presented 
in Figures U-lU through U-l8. In general the differences tend to 
be fractions of a count and are greater in the small areas 
(e.g. 2x2) where the calculations are more sensitive to Individual 
pixel differences. One except ■’on is in the high frequency region 
of scene B, band 5. This plot shows a higher standard deviation 
in the bi-resampled image. The nested square areas for which the 
mean and standard deviation are computed are centered on a sharp 
edge in the high frequency region. The greater standard deviation 
of the bl-resampled image is a manifestation of the overshoot of 
the cubic convolution algc. itfam when processing across an edge. 

An evaluation of this edge overshoot problem is presented in 
paragraph d) below. 

c) One-dimensional power spectra of scene A were computed 
for the three regions, both in the vertical and the horizontal 
directions. An 128-point Discrete Fourier Transform (DFT) was 
employed in the calculation of the spectrum of each line. Spectra 



of 50 contlfi^ouB lines were averaged by frequency in each of the 
regions. These averaged, or onoothed, spectra for the singly ?nd 
doubly resampled images in all three regions, plus the absolute 
difference of the singly and doubly resampled image spectra axe 
presented in Figures U-I 9 through U-36. In the computation of the 
spectra, the DC term was removed prior to the transform computation. 

The power spectrum response of the single and double 
resampling process is shown in Figure U-37. The frequency response 
of the bl-resampllng is simply the convolution of the single 
resampling response with itself. These curves suggest that, relative 
to the single resampling process, bl-resampllng boosts the lower 
frequencies and attenuates t.ie higher frequencies, with a crossover 
point of 0.29 cycles/sample, iln the Pigu' 3 **-19 through U-^b, this 
crossover frequency corresponds to I « 37.) For the high and medium 
frequency regions, and for both the horizontal and vertical spectra, 
the power spectra curves show this behavior. Crossover frequencies 


for these four cases are: 

Region Crossover 

high frequency, horizontal 37 

high frequency, verticed 1*3 

medium frequency, horizontal **1 

medium frequency, vert leal **5 


The deviations of the actual crossover points from the theoretical 
crossover points can be attributed to errors in estimating and 
smoothing the frequency response of the digital data, plus discretization 
errors in the radiometric count computed for each pixel in the 



and double resampll*' process. The low frequency regions 
do not exhibit this general behavior. This ip probably due to the 
computational errors, in the power spectrum estimation process, 
being larger than the actual energy present in the higher 
frequencies. 

d) The process of bi-resampling pri-'.arily degrades edges. 

This is due to the Suppression of the high frequency energy content 
in edges by the cubic convolution algorithm, as discussed above ir. 
c). This effect maniferts Itself as edge overshoot in the 
bl -resampled imagery. Several examples of this phenomenon are 
given In Figures «-38 through U-Ui. These figures compare the 
raaJometric values of the singly resampled imagery with those of 
the bl-resampled imagery along either a vertical or a horizontal 
line which crosses a sharp edge. Also included in these figures 
Is the "nearest neighbor" resampled raw data. (The rav data is not 
congruent to the r sampled data and is provided only as a guide to 
the geneial shape of the uncorrupted data. However, in Figure If -38, 
the raw data line deviates less than O.OU pixel from the 
resampled lines.) The first three figures cross an edge in the 
high freouency region of Hand County. The last figure crosses an 
edge containing a bright spot of a runway of Dulles Airport. 

These figures also show that there is no noticeable spreading of 
the edges due to the high frequency energy loss in the 
bi-resampled data. It should be r‘>emphasized h6i*e that the 
degradations, although measureable, are still slight and 
probably ignorable for most applications. 



This figure is a photograph contained 
in the pocket at the back of this report. 


Figure U-i. Scene E-1080-15192, Singly Resampled Image 
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Flfrure U-' 


This figure Is a photograph contained 
in the pocket at the heick of this report. 


. f^cene E-1080-15192, Doubly Resampled Image 
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This figure is a photograph contained 
in the pocket at the hack of this report. 


Figure U-3. Scene F-1080-15192* Difference Image 
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Pif?ure U-U(a) Scene E-1 080-15192 , Doubly Resampled Image, Histogram of High Frequency Region 
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Figure 4.14(B). Scene E-1080-15192, Medium Frequency Region, 
Radiometric Areal Statistics 
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Pljcure U-i6(a). Hand County, Band 5, High Frequency Region, 
Radiometric Areal Statistics 
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Figure U-17 (a). Hand County, Bamd 6, High Frequency Region, 
RadlOTietrlc Areal Statistics 
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Figure U-17(C). Hand County, Band 6, Low Frequency Region, 
Radiometric Areal Statistics 
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Figure 1»-18(A). Hand County, Band 7. High Frequency Re;-^lon, 
Radiometric Areal Statistics 
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Pifrure U-l8(B). Hand County, Band 7, Medium Frequency Region, 
Radioraet.'ic Areal Statistics 
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Flimre U-l8(C). Hand County, Band 7» Low Frequency Regi r.. 
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Figure U-26 Horizontal Spectrum, Low Frequency Region, 
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Figure U>27 Absolute Difference of Horisontal Spectra, 
Low Frequency Region, Scene E>1 080-13192 
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Figure h-29 Vertical Spectnm, High Frequency Region, 

Doubly Resaopled Image of Scene E>1060-1$192 
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Pl»rure U-36 


Absolute Difference of VerticeJ. Spectra, 
Low Frequency Region, Scene E-1 000-153 92 
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Figure 4-38. Edge Effects of Bi-resampling, Hand County, Band 5, 
High Frequency Region 


U-80 


i-j a a o r» 


raw data 

singly resampled data 



Vertical Data at Sample l60 


Figure U>39. Edge Effects of Bl-resampllng, Hand County, Band 
High Frequency Region 


U-81 



nMi:0i4H3OMO>S> 


raw data 


singly resaaqitled data 



Figure U-UO. Edge Effects of Bi>resanipling, Hand County, Band 5 , 
High Frequency Region 
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Vertical Data at Sample d2U 

Figure U-41. Edge Effects of Bi-resampling, Scene E-1080-1^192, 
Dulles Airport Region 
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r>ct1on 5. f^lanairicatlon Reeults 


*nii8 section discusses the results obtained In a comparison 
studv of the two resampled lAndsnt data sets of scene B ( Hand County, 

South Dakota ) over the LACTE Intensive Study Site. The imafre comparison 
was based on three analysis applications in IBM's Earth Resources Labortory 
(KRL): classification, image differencing, and bivariate histograming. 

The training and test fields used in classification and bivariate 
hi stogrammlng were checked against the ground truth data to identify the 
ground cover. The craining fields were grouped into the following classes: 
pasture, grass, oats, com, and spring wheat. In this discussion oft 
classification, channels 1 to U refer to bands U to 7 of the singly 
resampled version of scene B, and channels 5 to 8 refer to bands U to 7 
of the doubly resampled version. 

Separate classification processing using the maximum liklihood 
method was performed on each of the two images for the defined test and 
training fields. Two classification summaries were generated for each image, 
one with a 0% threshold and one with a 1% threshold for the training 
and test fields. These 0% and 1J» thresholds refer to the chi*squared 
test. Assuming the class statistics to be mitltivariate normal, the 1/K 
threshold, for example, means that the probability of a sample belonging 
to a class having a quadratic form value greater than chi -squared ( .01, U) 
is O.OI. 

The areal proportion estimate of each ground cover category 
Is computed as a function of the proportion of pixels assigned to the 
ground cover cilass. "Tie 1?5 threshold results Indicated that 3.U% of the 
pixels were unclassified for the singly resampled image and 5.5^ of the 
pixels were unclassified for the doubly resampled image. 
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l'h«> «r«al proportion estlmatas from the elaeslficatlon Bummarles show 
that the two Images correlate well with ground truth data. 


AREAL PROPORTION ESTIMATLG (?) 


IMAGE 

CORN 

OATS 

PAST/ 

SPRING 




GRASSieS 

WHEAT 

Single Resampled 

0% thresh. 

7.3? 

6.8? 

80.0? 

5.1? 

1? thresh. 

6.8? 

5.6? 

77.3? 

U.8? 

Double Resampled 

0? thresh. 

Y.5? 

6.7? 

79.8? 

5.8? 

if thresh. 

7.1? 

5.5? 

76.2? 

5.0? 

Ground Truth Map 

8.06? 

6M% 

81.23? 

U.2U? 


'''he results were derived from the classification of the ground truth 
test data. 

The bivariate histograms of the input data, Figxires 5-1 through 
5-H, show the variations in the values assigned to a given pixel by the 
resampling procedures. For example, in Figure 5-2 it can be seen that 
the pixels which were assigned a value of 29 for channel 2 ( band 5 in 
the singly resampled image ) were assigned values ranging from 25 to 32 
in channel 6 I band 5 in the doubly resampled image ) , and pixels which 
were assigned a value of 29 in channel C were assigned values ranging 
from 27 to 3^ in channel 2. There are 36 unique symbols {*,1,2, . . . ,9,A, 
B,C,...,Y,Z) used in the bivariate plot, with the later symbols in the 
aenuence Indicating a greater density than the earlier symbols. 

Inspection of the various bivariate plots shows that these high-density 
symbols are distributed along the diagonal of the plots, indicating that 
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the hulk of the pixeln in the singly and doubly resampled Images have 
the same value in both images. 

When the input channels are differenced and the difference 
images displayed ( Figures 5-5 to 5-8 ), no cleau* field boundaries or 
field structure is seen, which suggests that the pixels with different 
values in the two images are not related to the field structure. 

With the particular training fields defined in this study, 
the single resampled image classification results were closer to the 
ground truth data for three of the four classes. However, these differences 
in the classification results are slight, and cure within the acc\iracy 
limits of the classification algorithm. Therefore, the classification 
results Indicate no cleai preference for singly or doubly resampled data. 

The bivariate histogram for the two classification maps 
( Figure >-10 ), with channel one containing the singly resampled result 
ana channel two containing the doubly resamplea result, shows that 
approximately pixels, out of a total of 302U8, were classified 

the same in the two images. These numbers can be derived from the symbols 
on the diagonal. Therefore, approximately IIJS of the pixels were classified 
differently for the two Images. 

The classification difference map in Figure 5-9 shows the 
differing pixel ground cover assignments in the non-white areas. The 
white areas represent areas where the pixel values in the maps agreed. 

This map shows that the differences are scattered across the image, and 
do not appear to fall on field boundaries or within specific fields. 

In conclusion, it can be said that; 

® The positioning of the pixels I’hlch differ appears to be random 
and not associated ^d.th field boundaries for both the imagery 
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and the classification maps 

The classification results show no basis for concluding that one 
resampling procedure produces significantly less accurate 
classification results than the other. 
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Figure 5-3. Channel 3 and 7 Bivariate Histogram 
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Figure 5-6. Pixel Differences, Channels 2 and 6 
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Figure ;»-8. Pixel Differences, Channels U euid 8 
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APPKHDIX 


i. Xkiuations for project! ns' to the Lambert Conformal Conic Space 

equations mapping noints on the earth's surface expressed in 
geodetic latitude and longitude to the Lambert Conformal Conic (LCC) 
projection are presented here. 

A spheroidal earth model is assumed with semimajor and semiminor axes 

a = 6378165 

b = 6356783 


Basic parasieters of the ellipse are the eccentricity^ e, defined hy 

•7 


e a 


f/ 


la^ ♦ b^ 


and the principal radius of normal curvature, N, at latitude 
defined by 


“ / 5 — T" 

yi - eSin^^ 

Let and denote the two steuridca*d parallels of the projection. 

Let B ^2 )/2 and Aq denote the latitude euid longitude, 

respectively, of the origin in the LCC space. The mapping functions 
will be expressed in terms of the functions 

-iitan(Zj^/2) - i».tan(Z2/2) 
where Nj = N((^^), Z^ » Z«)j) 


tan( Z (^)/2 ) 


cot( fr/U ♦ ^/2 ) 


‘1 ♦ e sin<) 
.1 - e sln 4 > 


e /2 


A -1 



and 



N COB^i 

K ■ K(^, « — i- — 

^ ^ i (tan(2j^/2)) 

The maiiping equations are 

x(^,X) « r Bin$iX-^Q) 

y(^,X) « Tq - r coeJf(X-X^) 

where 

r » r(^,i]^,42) “ ^ tan^(Z(^)/2) 

and 

*“0 • *'< 4 ' o *^»^ 2 ^ 

Refer to Figure A-1 for the geometric interpretation of these 
relationships. 


2. Equations for mapping points frtns the I^C projection to geodetic 
latitude emd longitude in the northern latitudes. 

The equation for longitude, A , is straightforward. 

For latitude, it is necessary to solve 

/jC) « tan(fr/l» |/2) fl . e « f(^) 

Ir/ (iTVsiwq 

for This equation does not have a closed form solution. Therefore, 
an iterative procedure ( Newton-Raphson Method ) is used. 

Initially, ^ is set to a latitude corresponding to a spherical 
earth and then the Rewton>Raphson Method is applied to the parameter 



A-2 



where 


¥ 


f(^) • tan(ifA ♦ 


1 - e 3la^j^ 
1 ♦ e sin^ 


«/2 


Tn this method, the latitude is updated according to 

It is clear that finding a zero for the function P is equivalent to 
(^jj converging to 
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I'igure 3-i Feature Location, t'iuster Region, and 
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Figure J-2 Features from the Single Resanpled Imagery, 

Si'ene 1080 
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